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While not widely recognized, electronic correlations might play an important role in graphene. In- 
deed, Pauling's resonance valence bond (RVB) theory^ for the p7r-bonded planar organic molecules, 
of which graphene is the infinite extension, already established the importance of the nearest neigh- 
bor spin-singlet bond (SB) state in these materials. However, despite the recent growth of interest 
in graphene^ since its isolation,— there is still no quantitative estimate of the effects of Coulomb 
repulsion in either undoped or doped graphene. Here we use a tight-binding Bogoliubov-de Gennes 
(TB BdG) formalism to show that in unconventional d-wave contact graphene Josephson junctions 
the intrinsic SB correlations are strongly enhanced. We show on a striking effect of the SB correla- 
tions in both proximity effect and Josephson current as well as establishing a 1/(T — T c ) functional 
dependence for the superconducting decay length. Here T c is the superconducting transition tem- 
perature for the intrinsic SB correlations, which depends on both the effects of Coulomb repulsion 
and the doping level. We therefore propose that d-wave contact graphene Josephson junctions will 
provide a promising experimental system for the measurement of the effective strength of intrinsic 
SB correlations in graphene. 



Graphene is formed as a single layer of carbon atoms 
arranged in a honeycomb lattice with a peculiar band 
structure where the low energy excitations are massless 
Dirac quasiparticles which have a linear dispersion and 
a Fermi velocity vf ~ c/300. In comparison with quan- 
tum electrodynamics (QED), this low value for vf gives 
a large fine structure constant for graphene, a g ~ 2, 
and thus Coulomb interaction effects are expected to be 
strong. Recently this fact has been studied in terms of 
the possibility of Coulomb interactions driving undoped 
graphene insulating through correlations in the electron- 
hole channel^ but, in general, the effect of electronic 
interactions in graphene has mostly been ignored. The 
standpoint was, however, quite different in the early RVB 
treatments of the p7r-bonded planar organic molecules by 
Paulingi and others. Later, the RVB idea was also re- 
vived by Anderson^ as a mechanism for superconductiv- 
ity in the high-T c cuprate superconductors. BaskaranJ 
proposed in 2002 a Hamiltonian which phenomenologi- 
cally incorporates an effective JSi • Sj term (Si is the 
spin of atom i, and i and j are nearest neighbors), ul- 
timately derived from Coulomb interactions, as a model 
with which to estimate the SB correlations in graphene. 
A few years ago, we studied the Cooper pairing channel 
of this Hamiltonian and showed that a superconducting 
gap with d-wave symmetry develops below T c , where T c 
increases with the level of doping of the graphene. In the 
bulk, the (i-wave state takes the form d x 2_ y 2 + id xy and 
thus breaks time-reversal symmetry. Recently this state 
has been shown theoretically to survive on-site Coulomb 
repulsion^ as well as to emerge as an instability in a func- 
tional renormalization group flow, 9 results that both fur- 
ther strengthen the argument that superconducting cor- 
relations are important in graphene. 



There have been reports of granular superconduc- 
tivity in graphit e) 10 ' 11 ' 12 as well as a possible interfa- 
cial superconducting state between graphite and sulfur 
in the graphite-sulfur composites i 11 ' 13 ' 14 Both of these 
might be experimental realizations of the intrinsic SB 
mechanism,— but more experimental data is needed in 
order to establish such a conclusion. From an exper- 
imental point of view, the main obstacle to reach the 
SB superconducting regime in graphene is that a very 
heavy doping is needed for an experimentally measur- 
able T c . The most frequently used method of applying 
a gate voltage does not nearly induce enough carriers 
into the graphene whereas chemical doping often faces 
significant material science problems and can also pro- 
duce unwanted changes in other physical properties of 
the graphene. An exception to the latter possibly be- 
ing sulfur dopingii which offers an explanation to the 
experimental results for the graphite-sulfur composites 
mentioned above. Another viable approach would be to 
search for systems where the SB pairing correlations are 
naturally enhanced. One obvious candidate would be 
superconducting graphene. Induced superconductivity 
was recently demonstrated in graphene superconductor- 
normal metal-superconductor (SNS) Josephson junctions 
manufactured by depositing superconducting contacts on 
top of the graphene However, we recently showed that 
there is no enhancement of the SB correlations in conven- 
tional s-wave contact junctions^ The primary source 
for this behavior is the incompatibility between the d- 
wave symmetric SB state and the s-wave superconduct- 
ing contacts. It is therefore logical to expand our pre- 
vious study to include unconventional superconducting 
contacts where the symmetry of the contacts matches 
that of the SB superconducting state. Using the TB BdG 
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formalism we will report below a significant enhancement 
of the SB correlations in graphene d-wave contact SNS 
Josephson junctions. 

Method. — We start with the phenomenological Hamil- 
tonian for graphene proposed by Baskarani 
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Here /. and g. are the creation operators on the two 
inequivalent lattice sites of the honeycomb lattice, see 
Fig. [Hi, t = 2.5 eV is the hopping parameter, < i,j > in- 
dicates nearest neighbors, and \i is the effective chemical 
potential where fi = corresponds to the Dirac point. 
The SB correlations are modeled by the last term where 
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FIG. 1: (a) The graphene honeycomb lattice. Schematic ex- 
perimental setup (b) and model setup (c) for a d-wave contact 
graphene SNS Josephson junction. 

/jjj = — /i\.9j|)/V2 is the spin-singlet creation op- 

erator on the bond between atoms i and j. Note that 



J(Si • Sj 
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and therefore the last term in equation (jTJ) both effec- 
tively models the SB correlations and shows its Cooper 
pairing nature. The effective coupling constant J has 
been estimated to be as large as tZ Using Aj a (i) = 
i+ a ) as the mean-field order parameter, equa- 
tion UJ can solved with a double Bogoliubov-Valatin 
transformation. Here a denotes the nearest neighbor 
bond, see again Fig. [1^,, and we allow the order param- 
eter to be independent on the three bonds. The most 
favorable superconducting state in the bulk has A,/ a oc 
(1, e l27r / 3 , e* 4 ^/ 3 ) which corresponds to a d x i_ y i + id xy - 
wave symmetry in the Brillouin zone. As long as J > 
this state has a finite T c for any finite doping level and T c 
increases rapidly with either hole or electron doping^ 

In addition to equation (TT]), the influence of the super- 
conducting contacts, deposited on top of the graphene 
sheet as shown in Fig. [TJd, has to be taken into account. 
We do this by assuming that the contacts induce into the 
S regions of the graphene a) an effective paring potential 



and b) a high doping level. For d-wave contacts the sim- 
plest effective pairing potential is the same J-term as for 
the intrinsic SB correlations. We therefore chose to model 
the effect of the contacts with a strong J-term in the S re- 
gions, see Fig.[TJ:. The particular d-wave symmetry in the 
S regions is fixed to the d x 2_ y 2-w&ve on the zigzag inter- 
face which is achieved by setting Aj a oc (2,-1, — l) i 15 i 18 
We have checked that a 90° rotation of either interface 
or d-wave symmetry produces no notable change in the 
proximity effect. Note that the choice of d-wave sym- 
metry is still unrestricted in the N region, though with 
the symmetry fixed in the S regions and relatively short 
junctions, the same symmetry choice in N is now ener- 
getically favored. We should also note that we have set 
an artificially high value for J in S in order to get a 
short superconducting coherence length £ and therefore 
reduce the computational size of the problem. However, 
we believe our conclusions are still applicable for an ex- 
perimental setup with, e.g., a d-wave high-T c cuprate su- 
perconductor as the contact metal. The doping is set to 
be high in the S regions since the superconducting con- 
tacts act as large electronic reservoirs directly in contact 
with the graphene. In the N region we assume that the 
doping level is finite in order to produce a non-zero T c 
for the intrinsic SB correlations. In fact, the doping in 
the N region is kept moderately high in our calculations 
in order to deal with numerically accessible T c s for the 
SB superconductivity. But, since our results use only 
T c as the relevant SB parameter, this should not affect 
the generalization of our results to a more realistic lower 
doping regime. We assume that the effective chemical 
potential abruptly drops at the SN interfaces. While this 
is a simplification, a moderate leakage of charge from S to 
N would only raise the intrinsic T c in the interface region 
and thus enhance the SB pairing effect. 

The above framework can be solved self-consistently 
for the position-dependent Aj a (i) using the TB 
BdG formalism for ballistic graphene SNS Josephson 
junctionsj 17 ! 18 ' 19 We assume clean, smooth interfaces in 
order to reduce the computational effort but we do not 
expect a moderately disordered interface to significantly 
change our results^ It is the use of an effective pairing 
potential instead of an effective order parameter in the 
S regions that allows us to solve self-consistently for the 
order parameter. This leads to a self-consistent calcu- 
lation of the proximity effect depletion in S and leakage 
into N of the pair amplitude Fj a = Aj a /(v / 2J) and its 
norm F as seen in Fig. [5] This should be compared to 
the analytical Dirac-BdG (DBdG) formalism developed 
for graphene SNS Josephson junction s 20 ! 21 where the or- 
der parameter is assumed to be constant in the S regions. 
Not only is a self-consistent solution necessary for study- 
ing the coupling between the externally induced d-wave 
state and the SB pairing correlations but a recent com- 
parison between the two methods has also demonstrated 
the necessity for a self-consistent solution when uncon- 
ventional contacts are considered? 1 ^ 

Results. — The most straightforward physical quantity 
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to study in the TB BdG framework is the proximity ef- 
fect, or more specifically, the Cooper pair amplitude dis- 
tribution in the SNS junction. If a measurable signal of 
SB correlations were found in this quantity, there ought 
to be signs in other, more easily experimentally measur- 
able, quantities such as the Josephson current. Fig. [2] 
shows the norm of the pair amplitude, F, for a d-wave 
contact SNS junction. In black are the results when the 
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FIG. 2: Proximity effect in terms of the norm of the pair 
amplitude F normalized to the bulk value in the contacts Fo 
when SB correlations are ignored, i.e. J(N) = (black), and 
when SB correlations are present with J(N) = t (red), where 
t is the nearest neighbor hopping parameter. The results 
are given at temperatures T = 0.6T C (solid) and T = 17T C 
(dashed) where T c is the superconducting transition temper- 
ature for the intrinsic SB correlations. Vertical lines indicate 
the SN interfaces. The junction length is L — 30 unit cells 
to be compared with the superconducting coherence length 
£ = 5 unit cells in the contact regions. Inset shows a magni- 
fication over the N region. 

SB correlations are ignored at two different temperatures 
and, as seen, the results are essentially temperature in- 
dependent in this range. When SB correlations are in- 
cluded (red) the pair amplitude in N increases. The pair 
amplitude inside N is in fact larger by almost an order of 
a magnitude than it would have been in a bulk sample 
with the same strength SB pairing. This means that in 
d-wave contact SNS junctions the intrinsic SB pairing is 
heavily benefitting from the proximity to the supercon- 
ducting contacts. This is despite the large Fermi level 
mis-match (FLM) at the SN interfaces, due to different 
doping levels in S and N, which acts as an effective in- 
terface barrier»iS However, since the extrinsic supercon- 
ducting state and the intrinsic SB pairing here by design 
have the same symmetry this is probably not a too sur- 
prising result. More interesting is the temperature scale 
over which this enhanced SB pairing takes place. The 
solid curve is for T — 0.6T C , i.e. when the doping level is 
sufficient for intrinsic SB superconductivity to be present, 
and the dashed curve is recorded when T = 17T C . We 
thus see a strongly enhanced proximity effect in the SNS 
junction even far above T c . 

One very closely related, but easier quantifiable, prop- 
erty to the proximity effect is the superconducting decay 



length, or normal state coherence length £„, which is the 
length scale over which the pairing amplitude (exponen- 
tially) decay inside N: 



F oc e 



(3) 



Fig. [3] shows the inverse temperature dependence of £ n 
for several different strengths of the intrinsic SB pairing. 
In black circles are the results when SB correlations are 
ignored whereas colored crosses show the results when 
they are present with the vertical lines indicating the 
corresponding intrinsic T c . The red and magenta curves 
correspond to the same T c but achieved through differ- 
ent combinations of SB correlation strength J and dop- 
ing level, which demonstrates the dependence on T c only. 
From the figure we clearly see that to a good approxima- 
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FIG. 3: Superconducting decay length £ n as a function of 
scaled inverse temperature t/(ksT) for J(N) = (black, cir- 
cles) and J(N) > (colors, crosses) with increasing transition 
temperatures T c achieved by increasing J(N) (yellow, cyan, 
blue, red, green). Increasing the doping level in N instead 
(magenta) produces the same result. Vertical lines indicate 
the corresponding T c s. The inset shows the linear behavior of 
£„ vs. t/(k B {T - T c )) at moderately to large (T - T c ). 

tion £„ cx 1/T in the uncorrelated case. For any finite 
SB pairing inside N, £ n will diverge at the correspond- 
ing intrinsic T c leading to a giant proximity effect for 
temperatures close to T c . The same diverging behavior 
has recently also been predicted for the square lattice 
for both s-wave and ci-wave contact SNS junctions when 
the N region exhibits weak superconductivity of the same 
symmetry^ although this study was conducted with no 
FLM at the SN interfaces. The increase in proximity ef- 
fect far above T c documented in Fig. [2j is here reflected 
in an increased decay length even at elevated tempera- 
tures. The inset in Fig. [3] shows that for temperatures 
moderately to high above T c we can extract the following 
functional form for the superconducting decay length: 
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FIG. 4: Critical current I c in units of Io = eW/ (H£) when 
vf is set to 1 for J(N) = (circles) and J(N) = t (crosses). 
Different colors represent increasing doping levels (black, red, 
green) in N. The corresponding T c s are marked with vertical 
lines. The junction length is L — 30 unit cells to be compared 
to £ = 5 unit cells. 

The limited accuracy due to finite model systems makes 
our numerical data deviate from this relationship at tem- 
peratures close to T c , but we expect equation (0| to be 
valid even as T — > T c . From this we can conclude that 
any evidence of a divergence in £„ at a non-zero tempera- 
ture with a 1/(T—T C ) functional form in a d-wave contact 
graphene SNS junction would be a clear measurement of 
the strength of the intrinsic SB correlations in graphene. 

It would also be useful to characterize the effect of SB 
correlations on the most accessible property of a SNS 



junction, namely the Josephson current. Fig. [4] shows 
the critical current as a function of temperature when SB 
correlations are ignored (circles) and present (crosses) for 
three different doping levels in N. The corresponding T c s 
for the SB pairing are again marked with vertical lines. 
As seen, the Josephson current increases with decreasing 
temperatures in all cases. This is expected since £ n is 
inversely proportional to the temperature. There is also 
a pronounced increase in the current when SB correla- 
tions are present in the junction, and, most importantly, 
the current is enhanced far above T c . For systems with 
stronger SB correlations, the current increases more com- 
pared to the uncorrelated case at and below T c . This 
effect is not due to the decay length per see, since £„ di- 
verges at T c for all strengths of the SB correlations, but 
related to the strength of the intrinsic SB superconduct- 
ing state in N. 

Discussion. — The results in Figs. [2E1 clearly demon- 
strate that the intrinsic SB correlations in graphene are 
greatly enhanced in a graphene d-wave contact SNS 
Josephson junction even far above their intrinsic T c . 
Since the intrinsic SB correlations at the currently ex- 
perimentally fairly limited doping levels of graphene are 
expected to be rather weak, this enhancement even far 
above T c is crucial for any experimental detection and a 
finding of a finite T c in a graphene d-wave contact SNS 
Josephson junction will directly establish the strength of 
the effective SB coupling J as a function of the doping 
level of the graphene. 
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